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Abstract 

When described in a grand canonical ensemble, a finite Coulomb system exhibits 
charge fluctuations. These fluctuations are studied in the case of a classical (i.e. 
non-quantum) system with no macroscopic average charge. Assuming the validity 
of macroscopic electrostatics gives, on a three-dimensional finite large conductor of 
volume V, a mean square charge < > which goes as V^^^. More generally, in 
a short-circuited capacitor of capacitance C, made of two conductors, the mean 
square charge on one conductor is < >= TC, where T is the temperature and 
C the capacitance of the capacitor. The case of only one conductor in a grand 
canonical ensemble is obtained by removing the other conductor to infinity. The 
general formula is checked in the weak-coupling (Debye-Hiickel) limit for a spherical 
capacitor. For two-dimensional Coulomb systems (with logarithmic interactions), 
there are exactly solvable models which reveal that, in some cases, macroscopic 
electrostatics is not applicable even for large conductors. This is when the charge 
fluctuations involve only a small number of particles. The mean square charge on 
one two-dimensional system alone, in the grand canonical ensemble, is expected to 
be, at most, one squared elementary charge. 
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1 INTRODUCTION 



This paper is dedicated to Michael Fisher on the occasion of his 70th birthday. Part of 
it deals with the Debye-Hiickel theory of Coulomb systems, to which Michael brought 
elaborate refinements, in particular insisting on the importance of taking hard cores into 
account. I apologize for using here only the simple point-particle version. 

I consider the classical (i.e. non-quantum) equilibrium statistical mechanics of Coulomb 
systems: systems of charged particles interacting through the Coulomb law (plus perhaps 
short-ranged forces), such as plasmas or electrolytes. Such a system, when described by 
a grand canonical ensemble, is expected to exhibit charge fluctuations. The aim of the 
present paper is to study these fluctuations, for a finite but macroscopic system. 

The grand canonical ensemble is often introduced by considering a system in contact 
with a surrounding infinite system (the reservoir), with possible exchanges of energy and 
particles between the system and the reservoir. In this approach, it is assumed that the 
energy of interaction between the system and the reservoir is negligible. This is indeed the 
case for a finite but macroscopic system, when the interparticle forces are short-ranged 
(then, the system energy goes as its volume, while the interaction energy goes only as 
its surface area). However, this approach has to be modified when there are long-ranged 
forces such as Coulomb ones. Then, for the interaction energy between the system and 
the reservoir be disregarded, is is necessary to assume that the reservoir is infinitely far 
away from the system under consideration. This is how the grand canonical ensemble will 
be defined in the following. 

Studying the charge fluctuations in a given large subvolume A of an infinite Coulomb 
system (i.e. assuming that the reservoir is in contact with the subvolume) is a different 
problem, which has already been studied and solved[^0. It was found that the mean 
square charge (Q^) in A behaves as its surface area S (not its volume V). In the presently 
studied geometry (infinitely remote reservoir), it will be argued that {Q^) is even weaker, 
behaving as V^^^, for a 3-dimensional system. This behavior had been conjectured by 
Lieb and Lebowitz0, but they could not prove it rigorously, i.e. by a purely microscopic 
argument. Here, on the contrary, our starting point will be the validity of macroscopic 
electrostatics of conductors, assumed without proof. 

In general, in the grand canonical ensemble, the state of a system, made of s species 
of particles, depends on s chemical potentials. However, in the case of a Coulomb system, 
in the thermodynamic limit, the state depends only on s — 1 chemical potentials and is 
the same as in a grand canonical ensemble restricted to neutral configurations [Q]. In the 
case of a macroscopic but finite Coulomb system, there may be a non-zero average charge 
(Q) and one more chemical potential is needed for controlling it. Here, we only consider 
the simple case in which there is only one reservoir, which is an infinite Coulomb system 
of the same nature as the system under consideration, and (Q) is expected to vanish if 
macroscopic electrostatics holds. 

We shall first assume that our system and this reservoir are at some distance of each 
other. Since they can freely exchange particles, the system and the reservoir can be 
considered as two conductors forming a short-circuited capacitor. In Section 2, a general 
simple expression for the mean square charge of the system in that configuration will 
be derived; the case of the system in the grand canonical ensemble will be obtained by 
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removing the reservoir to infinity. In Section 3, these general results will be checked by a 
microscopic calculation in the weak-coupling (Debye-Hiickel) limit. 

Coulomb systems can be mimicked in a two-dimensional world, in which the Coulomb 
potential 1/r must be replaced by the two-dimensional solution of the Poisson equation 
— Inr. Working in two dimensions has the advantage that exactly solvable models for 
the statistical mechanics of Coulomb systems are available. However, in two dimensions, 
some specific subtle points arise and deserve a separate discussion. In Section 4, two- 
dimensional models will be considered, and cases when macroscopic electrostatics cannot 
be used will be exhibited. 

Rather than starting with a capacitor made of two conductors and removing one of 
them to infinity, one might want to study directly the case of one finite system in the 
grand canonical ensemble. How to formulate the Debye-Hiickel theory in a finite system 
is discussed in Appendix B. 



2 CHARGE FLUCTUATIONS IN A CAPACITOR 

Let us consider a capacitor, made of two macroscopic conductors A and B separated 
by vacuum, with B surrounding A. Let C be the capacitance. Let Q be the charge on 
conductor A, —Q the charge on conductor B. When the capacitor is short-circuited, the 
average charge (iQ) on each conductor vanishes. If the conductors are describable by 
classical (non-quantum) statistical mechanics, we claim that the mean square charge on 
each conductor is 

{Q^)=TC (2.1) 

where T is the temperature (in units of energy). This relation is a special case of Nyquist's 
formulaic which gives the electrical current fluctuations in a linear electric circuit. The 
derivation of Nyquist's formula relies on the fluctuation-dissipation theorem in the theory 
of linear response. In the present case, a more direct derivation of eq.( |2.1| ) can be made 
by using the simpler classical static linear response theory, as follows. 

Let us introduce some given external charges which create an infinitesimal electric 
potential difference between the two conductors: let the potential on conductor A minus 
the potential on conductor B be 5(j) (for instance, in the case of a spherical capacitor, 
we introduce between the conductors two spherical concentric layers carrying opposite 
uniformly distributed charges; then, these layers create a potential which is zero on the 
outer conductor B and has some constant value S(j) on the inner conductor A). The 
corresponding change in the Hamiltonian is H' = S(j)Q, where Q is the charge on conductor 
A. The short-circuited capacitor will respond by tranfering a charge 6Q = —C 6(j) from 
conductor B to conductor A, in such a way that the total potential difference (external 
plus induced) between the conductors vanishes. Now, linear response theory says 

-C6(f) = 6Q = -(3{{H'Q) - {H') (Q)) = -fi 5<P{Q^) (2.2) 

where is the inverse temperature 1/T and the averages (. . .) are taken in the absence of 
the perturbation H' (we have used {Q) = 0). Eq.(|2.2|) proves eq.(0|).lt should be noted 



that C will be the usual capacitance for a given geometry only provided that macroscopic 
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electrostatics is applicable. A necessary condition is that the sizes of the conductors and 
the separation between them be large compared to the microscopic scale. 

The mean square charge of a finite macroscopic Coulomb system in the grand canonical 
ensemble is obtained by sending conductor B to infinity. Then the mean square charge 
on the macroscopic body A is given by ( |2.1| ), where now C is the capacitance of the 
macroscopic body A alone. For a 3-dimensional system, this capacitance goes as V^^^ 
where V is the volume, as stated in the Introduction. 

In the special simple case of a spherical capacitor, made of an inner conductor of radius 
-Ri and an outer conductor of radius R2, the capacitance is 

It becomes i?i if i?2 C)0, giving the mean square charge in a macroscopic but finite 
spherical Coulomb system of radius R, in the grand-canonical ensemble, 

(Q2) = TR (2.4) 

in agreement with previous findings about the charge correlations. |^ 

Another limit of interest is when Ri becomes large compared to W = R2 — Ri which 
keeps a fixed (macroscopic) value. One then obtains a plane capacitor with plate areas 
S = inRf and plate separation W, and indeed ( p.3| ) becomes C = S/i^AnW), giving for 
the mean square charge per unit area on one plate 

Un, M = ^ ,2.5) 

5-00 S 47lW ^ ' 



3 WEAK-COUPLING LIMIT 

We consider the simple geometry of a spherical capacitor made with a classical 
Coulomb fluid: There are two concentric spheres, centered at the origin, of radii -Ri 
and i?2 (-R2 > Ri)- The shell between the spheres is empty. The ball of radius Ri is 
filled with the fluid, as well as the whole space outside the sphere of radius i?2- The 
short-circuiting of the capacitor is described by assuming that the two filled regions are 
allowed to freely exchange particles. 

The weak-coupling limit is a high-temperature one which is expected to be correctly 
described by the Debye-Hiickel theory. The Coulomb fluid is made of several species of 
particles of number densities Ua and charges qa- For the system to be stable, in addition 
to the Coulomb forces, there should be some short-ranged repulsive forces, but the weak- 
coupling limit can also be viewed as a low-density one in which these short-ranged forces 
can be neglected. Strictly speaking, the number densities are position-dependent near the 
fluid boundaries. However, taking this into account would only give corrections of higher 
order and therefore we consider the densities as constants. The Debye wave number is 
defined as k = (47r/5 '^a^a)^'^^- Let p(r) be the microscopic charge density at r. We 
shall need the charge correlation function 

(P(r)p(r )) = E naqlnmlKiv, r') + ^ n,ql5{v - r') (3.1) 

ab a 
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where both r and r' are in a filled region. K{r, r') is the solution of the partial differential 
equation 

[A - K'^{r)]K{r, r') = A7ipS{r - r') (3.2) 

where the source point r' is assumed to be in a filled region, while r can be anywhere: 
K^(r) = if r is in a filled region and /t^(r) = if r is in the empty region. K and its 
normal derivative must be continuous at the boundaries r = Ri and r = R2, and K 
as r — > cx). 

In the present Debye-Hiickel scheme, the average charge on the sphere of radius Ri 
vanishes, while the mean square charge on that sphere is 

iQ') = I dr / dr'(p(r)p(r')) (3.3) 

Jr<Ri Jr'<Ri 

The solution K of ( p^ ) is studied in Appendix A. When used in ( p.l| ) and ( |3.3| ) it gives 

^^g2\ _ [1 + kR2][hRi cosh(/t-Ri) - sinh(ft:i?i)] 

k[{1 + KR2 — K,Ri) cosh(fi;_Ri) + sinh(/ti?i)] 



If both Ri and R2 are macroscopic, i.e. if kRi, KR2 ^ 1, ( p^ becomes 

The mean square charge in a large spherical subdomain of an infinite Coulomb system 
is retrieved by taking Ri = R2 = R 'vn (p.5|) which becomes 



/3(Q2) = ^t,S (3.6) 
where here S = AirR"^ is the sphere area, in agreement with the general formula[|l| 



-i|drr(p(0)p(r)) (3.7) 



S 

where (p(O)p(r)) is the infinite-system charge correlation function, obtained by using in 
( p.2[ ) the infinite-system function i^'(0,r) = — /5exp(— /s:r)/r. 

If, on the contrary, R2 — Ri is also macrosopic, i.e if we now assume k{R2 — Ri) ^ 1 as 
well as kRi ^ 1 in (|3.5D, we do check the expected spherical capacitor charge fiuctuation 



!T]) with capacitance (2.3), as well as the grand-canonical fiuctuation ( |2.4| ). The limit 
( ^.5[ ) of a plane capacitor can also be taken. 

The plane capacitor can also be studied directly. There are two parallel planes a; = 
and X = W, thus separated by a distance W. The slab between the plates < x < W 
is empty, while the Coulomb fiuid fills the two semiinfinite regions x < and x > W 
outside the slab. In this geometry, it is possible to solve ( p.2|) for a function K which now 
is a function of x, x', and the component y of r — r' along the plates (actually, one rather 
computes the Fourier transform of K with respect to y). Then, one obtains the mean 
square charge per unit area on one plate of (infinite) area S as 



lim ^ = [dy [ dx [ dx'(p(a;,y)p(x',0)) (3.8) 

5^00 J Jx<0 Jx'<0 
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The result is 

s^oo S A'k{2 + kW) ^ ' 



If = 0, one retrieves (|3.6|) . If W is macroscopic, i.e. kW ^ 1, one retrieves 

The above considerations also apply to the case of a one-component plasma (also called 
jellium), i.e. a system made of mobile point charges of one species, with number density 
n and charge g, in a uniform neutralizing background. The Debye wave number reduces 
to K = (47r/?ng^)^/^. However, for classical jellium, the grand canonical partition function 
is a convergent series only if the background is kept fixed while one sums over the particle 
number 0. The same kind of prescription should be used here: regions A and B are 
assumed allowed to exchange particles, but the backgrounds keep a fixed charge density. 



4 TWO-DIMENSIONAL COULOMB SYSTEMS 

4.1 General Properties 

In the two-dimensional systems discussed in this Section, the electric potential created 
at r by a unit charge at the origin is — ln(r/L) where L is some fixed length. This choice 
of a two-dimensional solution of the Poisson equation often makes these systems good toy 
models for mimicking three-dimensional systems with the usual 1/r potential. One of the 
advantages of working in two dimensions is the existence of exactly solvable models. For 
avoiding any confusion, it should be stressed that these toy models do not describe "real" 
charged particles such as electrons, which, even when confined in a surface, still interact 
by the 1/r law. 

In the simple case of a circular capacitor, made of an inner circular conductor of 
radius i?i and an outer circular conductor of radius -R2, macroscopic two-dimensional 
electrostatics says that the capacitance is 



C goes to as i?2 — * OO.0 More generally, the macroscopic capacitance of one finite 
conductor of any shape vanishes. This is because bringing from infinity an additional 
charge q onto a conductor of characteristic size i?, already carrying a charge Q, would 
cost an energy of order qQ dr/r which is infinite. 

It will however be shown below that, for the present problem of charge fluctuations, 
the mimicking of three-dimensional systems which might be expected at first glance does 
not always occur, because, in some cases, although all relevant lengths are macroscopic, 
nevertheless macroscopic electrostatics cannot be used. In the case of a circular capacitor, 
(Q^) may not be given by using (|4.1| ) in (|2.1|) . 



^In a previous paper ||^, another definition of the capacitance of a disk was used. The present one 
(the limit of the capacitance of a circular capacitor when the outer conductor recedes to infinity) is more 
appropriate here. 
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4.2 Weak- Coupling Limit 

The two-dimensional case is very similar to the three-dimensional one. We now 
consider the simple geometry of a circular capacitor replacing the concentric spheres by 
concentric circles: The regions r < Ri and r > R2 are occupied by the Coulomb fluid, 
while the region Ri < r < R2 is empty. The Debye-Hiickel theory is again used, with Att 
replaced by 27r in the definition of the Debye wave number k and in the r.h.s. of ( |3.2| ). 
The detail of the calculation is given in Appendix A. The result is 

h(KRi) R2 Ki(kR2) '^-^l Ri 

where J/ and Ki are modified Bessel functions, while (Q) = 0. If both Ri and R2 are 
macroscopic, i.e. if kRi, KR2 ^ 1, (|4.2|) becomes 

f3{Q^) = (4.3) 

The mean square charge in a large circular subdomain of an infinite system is retrieved 
by taking Ri = R2 = R in ( |4.3| ) which becomes 

m') = (4.4) 

where S = 2ttR is the subdomain perimeter, in agreement with the two-dimensional 
analog! of Q: 

m 1 



j drr(p(0)p(r)) (4.5) 



S vr 

where {p{0)p{r)) is the infinite-system charge correlation function, obtained by using in 
( p.2| ) the infinite-system function K{0,r) = —PKo^Kr). 

If, on the contrary, R2 — Ri is also macroscopic, i.e. if k{R2 — Ri) ^ 1 as well as 
uRi ^ 1, ( [4.3| ) becomes 

m') = (4.6) 

in agreement with (|2.1| ) and ( [4.1| ). In the limit R2 — > 00, there is no charge fluctuation. 

4.3 Two- Component Plasma at F = 2 

The two-dimensional two-component plasma is made of two species of point particles 
of opposite charges ±g, interacting through the pair potential ±g^ln(r/L). The dimen- 
sionless couphng constant is F = Pq"^. The system is stable against collapse of pairs of 
oppositely charged particles for F < 2. Many exact results are now available for that 
model in its whole stability range [|^ 11 1. However, the correlation functions are fully 



known only at F = 2, for a variety of plane p^-[T7| (or even curved) geometries (although, 
for a given fugacity, the density starts to diverge at F = 2, the truncated many-body 
distribution functions remain finite). 
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By an extension of previous calculations, the correlation functions at F = 2, in the 
circular capacitor geometry described in Section 4.2 are obtained in Appendix C. The 
fugacity z (which is the same for both species) appears through a parameter m = 2TiLz 
which has the dimension of an inverse length and defines a microscopic scale (the bulk 
correlation length is (2m) ~^). Using these correlation functions in ( p.3| ) gives 

' \lIi{X)Ki{aX){a' - a~') - Xk{X)K[{aX)a'+^ + XI[{X)Ki{aX)a-^]' 

(4.7) 

where X = mRi and a = i?2/-Ri- By charge symmetry, (Q) = 0. 

The mean square charge in a large circular subdomain of radius R of an infinite 
system is retrieved by taking X = mR and a = 1 in ( |4.7| ) and using for the modified 



Bessel functions Ii and Ki the leading terms of their uniform asymptotic expansions |18 
which are appropriate in the present case of a large argument X and an index / which 
may also be large. It is found that the sum on I can be replaced by an integral. The 
result is 

{Q') = g'jS (4.8) 
where S = 2ttR is the subdomain perimeter, in agreement with the general formula ([4.5|), 



where the infinite-system charge correlation function ||13|,|14| (p(O)p(r)) is (up to its here 
irrelevant 6 term) —2q'^[m'^/{27r)]'^[KQ{mr) + Kf{mr)]. 

If, on the contrary, not only Ri but also R2 — Ri are macroscopic, i.e. X ^ 1 with 
a > 1, because of the a', a~\ a'"*"^ terms in the denominator of ( [4. 71) only small values of 
\l\ contribute to the sum, and the modified Bessel functions can be simply replaced by 



the leading terms of their ordinary asymptotic expansions [p^ , valid for a large argument 
X or aX and a given index /. In the limit X = mRi ^ 00 at fixed a = R2/R1, with 
= 2 being taken into account, (|4.71 ) becomes 



00 



i=-oo (^a'+2 + a-('+2) j 

where a = i?2/-Ri- It is clear that, when a has a finite value, the sum in ( [4.9| ) cannot be 



replaced by an integral [which would reproduce (^l6|) ]. 

( [4.9| ) is, at first sight, a very surprising result. It does not reproduce the value ( [4.6| ) 
expected on the basis of macroscopic electrostatics. On second thought, one sees that 
the l.h.s. of ( [4. 91 ) can be rewritten as 2{Q'^)/q'^ since here F = /?g^ = 2, and therefore 
{Q'^)/q^ is of order unity, which means that the fluctuation involves only a small number 
of particles. Thus, in spite of the fact that the relevant lengths Ri and R2 — Ri are 
macroscopic, the number of involved particles is not, and (|4.6|) based on macroscopic 
electrostatics should not be expected to hold at F = 2, and more generally at any finite 
temperature. However, in the weak-coupling (i.e. high-temperature) limit considered in 
Section 4.2, Pq"^ — > 0, as given by ( |4.6| ) becomes large, and this result is consistent 

with macroscopic electrostatics, as it should. 

It should be noted that, like (|4.6| ), (^]9|) indicates that there is no charge fluctuation 
in the limit R2 00, i.e. for one macroscopic disk in the grand canonical ensemble. Only 
the charge Q = contributes to the grand canonical distribution. 
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Another limit of interest is when Ri ^ oo for a fixed value ofW = R2 — Ri- One then 
obtains a two-dimensional "plane" capacitor. When this limit is approached, Ina ~ W/Ri 
is small, and the sum in ([4.9|) can be replaced by the integral 



1 dl 27ri?i 



This is the result expected on the basis of macroscopic electrostatics, the two-dimensional 



analog of (|2.5| ) with now the plate area replaced by a plate length 2ttRi and the capacitance 
the two-dimensional one C = 2ttRi/{27iW). Since {Q^)/q'^ now becomes large as Ri 00, 
macroscopic electrostatics should indeed hold. A direct derivation for a plane capacitor 
is feasible. 



4.4 One- Component Plasma at F = 2 

The two-dimensional one-component plasma is made of one species of point particles of 
charge q, interacting through the pair potential —q^ ln(r/L), in a neutralizing background 
of fixed charge density —qn. Far from the boundaries, the particle number density is n. 
The dimensionless coupling constant again is F = Up to now, the system is exactly 
solvable ||19|| only at F = 2, in which case the correlation functions are known for a large 
variety of plane pO|-|23| (or even curved) geometries. 

By an extention of previous calculations, the charge density and the correlation func- 
tion at F = 2, in the circular capacitor geometry, are obtained in Appendix D. The 
background is fixed, while the two regions can freely exchange particles. There is no 
charge symmetry and the average charge (Q) on the inner disk does not automatically 
vanish. Using in (|3.3|) (modified for taking (Q) into account), the correlation function of 
Appendix D gives the charge fluctuation on the inner disk. We introduce the notations 
Yi = irnRf and Y2 = TmR\ and the incomplete gamma functions [p^ 



and 



7(/+ l,ri) = / 'dte-V (4.11) 
Jo 

T{l + l + Y2-Yi,Y2) = / dt e'V+^^~^^ (4.12) 

JY2 

The average charge on the inner disk is found to be 

7(/ + 1, Yi) + D{L + 1, Yi, Y2) 

where 

D{1 + 1, Fi, Y2) = exp(Fi InFi - Yi - Y2 InFa + Y2)T{1 + 1 + ^2-^1, Y2) (4.14) 
The charge fluctuation on the inner disk is found to be 

w->-w>'^.-£ ,:'!r!;!^'^';/;^^.,. (4.15) 



-[7(z + i,ri) + D(/ + 1,^1,^2 



2r' 
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The integrands in ([4 .111 ) and ( |4.12|) have a maximum (a sharp one when / is large) at t = Z 
and t = I + Y2 — Yi, respectively. Therefore, 7(/ + 1, Yi) becomes small when / > Yi, while 
r(/ + 1 + 1^2 ~ ^1, Y2) and thus D{1 + l,Yi, Y2) become small when I < Yi. The summand 
in (|4.15| ) has a maximum near I = Yi. 

The results for a large circular subdomain of radius R of an infinite system are retrieved 
by taking = = R, i.e 1^1 = ^2 = TmR\ Then (|41^ ) gives (Q) = 0. Using for the 



incomplete gamma functions ( [4.11| ) and ( [4.12[ ) the Tricomi asymptotic representation 



7(/ + 1, r) ~ r(/ + l)[^ + ^-'/'Erf(^^^)] (4.16) 

where F is the complete gamma function and Erf is the error function (this representation 
is appropriate when Y is large and / is close to Y), it is found that the sum on / in (|4.15|) 
can be replaced by an integral. The result is 

iQ') = e'^S (4.17) 
Ztx 

where 5* = 2ttR is the subdomain perimeter, again in agreement with the general formula 
( [4.5[ ), where the infinite-system charge correlation function pO[] (p(O)p(r)) now is (up to 
its here irrelevant 5 term) — g^n^ exp(— vrnr^). 

We now turn to the case when not only Ri but also R2 — Ri are macroscopic, i.e. 
Yi ^ 1 with a = R2/R1 = (^2/^1)^^^ > 1- The average charge (|4.13 ) does not seem to 



have a simple expression. For investigating the charge fluctuation ( [4. 15] ), it is convenient 



to rewrite the summand in it as [(1^/7)^/^ + (7/D)^/^] ^. As it will be seen a posteriori, 
now the relevant values of / — Yi in the sum ([4.15|) are only of the order of a few units 
(rather than of the order of Y^^"^) and therefore the Erf term in the Tricomi representation 
( [4.16| ) can be omitted and the incomplete gamma functions just replaced by half the cor- 
responding complete one. Furthermore, one can use for these complete gamma functions 
the Stirling asymptotic expansions F(l + /) ~ exp[/ In / — Z + (1/2) ln(27r/)] and the similar 
one for F(/ + 1 + 1^2 ~ ^i)- Neglecting terms of order l/Yi or smaller, for a given value of 
/ — Yi, in ln(D/7), one finds 

Finally, one can shift the summation index in ( [4. 151) , replacing / — Yi by / — Fi , where Yi 
is the non-integer part of Yi such that < Y^ < 1, and extend the summation to —00 in 
the present large- Yi limit. Thus, j3q^ = 2 being taken into account, 

i=-oo (a'"^i+2 + a-('-^i+2) j 

where a = R2/ Ri- The form of ( [4.19| ) justifies a posteriori our above statement that the 
relevant values of / — Yi in ( 4.15| ) are only of the order of a few units. 



Here too, for a finite value of a, the sum ( [4.19| ) cannot be replaced by an integral and 



does not reproduce the value (^^) expected on the basis of macroscopic electrostatics. 
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The reason is the same as in the case of the two-component plasma: the fluctuation 
involves only a small number of particles. If Yi = 0, i.e if the background charge —qnirR^ 
on the inner disk is an integer number of elementary charges — g, the fluctuation ( |4.19|) 
for the one-component plasma is the same as the fluctuation (|4.9|) for the two-component 
plasma; we have no explanation to offer. 

The case of one disk alone, in the grand canonical ensemble (with a fixed background), 
is obtained by taking the limit R2 ~^ oo. Now the average charge on the disk has a simple 
form. For obtaining it, it is convenient to rewrite the summand in ( 4.13 ) as [1 + [D j^)]^^, 



where now, from (^^, D/^i = if / < Yi - (1/2), D/^i = 1 if / = - (1/2), and 

fefore, since the sumn 

(g) = -gFi ifFi< 



D/7 = cxD if / > Yi — (1/2). Therefore, since the summation index I is an integer 

1 



2 

(Q)=0 ifFi = i 

(Q) = g(l-Fi) if?i>^ (4.20) 



For the behavior of the fluctuation ( |4.19| ), in the limit a —>■ 00, two cases have to be 
distinguished. 

{Q')-{Qr = ifFi^l 

{Q')-{Q)' = j if>l = ^ (4.21) 

These results have a simple interpretation: —qYi = —qimRl is the negative background 
charge. Taking into account a (necessarily integer) number of positive particles gives a 
total charge (background plus particles) Q. If Yi ^ 1/2, only one value of Q contributes 
to the grand canonical ensemble: the one which corresponds to the smallest possible value 
of \Q\. If however Yi = 1/2, this smallest possible value is \Q\ = {q/2) which corresponds 
to two possibilities Q = ±(g/2) with equal probabilities. 

The two-dimensional "plane capacitor" limit (-Ri 00 for a fixed value of W = 
R2—R1) is the same as in the case of a two-component plasma. As this limit is approached, 
the sum in ( [4.19| ) can be replaced by an integral, and {Q) = because of the geometrical 
symmetry between the two flat plates. Again, the result is (|4.10| ), in agreement with 
macroscopic electrostatics. 



5 SUMMARY AND CONCLUSION 

For studying the charge fluctuations on a macroscopic but finite classical Coulomb 
system in the grand canonical ensemble, i.e. when the system is allowed to exchange 
particles with a reservoir, in a first step we have considered a capacitor: one electrode is 
the finite system under consideration, the other electrode surrounds the first one at some 
distance of it and extends to infinity. Both electrodes are assumed to be made of the same 
Coulomb fluid. The capacitor is short-circuited, which means that the two electrodes can 
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freely exchange particles. When the external electrode recedes to infinity, the internal one 
becomes one Coulomb system in a grand canonical ensemble. For actual calculations, the 
simple geometry of a spherical capacitor has been chosen. 

For a three-dimensional system, there is no surprise. A short-circuited capacitor of 
capacitance C can be considered as an electric oscillator, and it is rather natural to 
state that its average energy {Q'^)/2C is (1/2)T: this gives ( |2.1| ). This general formula 
is supported by the derivation of Section 2, on the basis of linear response theory and 
macroscopic electrostatics. It has been checked in Section 3, in the Debye-Hiickel theory. 

Two-dimensional Coulomb systems (with a logarithmic interaction) are more tricky. 
Since some of them are exactly solvable, it was tempting to test on them the general 
formula ( |2.1|) for the charge fiuctuations. It has been a surprise for the author that this 
general formula is not valid for a circular capacitor at some finite temperature. On second 
thought, one realizes that the charge fiuctuations involve only a small number of particles, 
and therefore one should not expect the validity of a macroscopic formula. 

In the limiting case of one disk alone in a grand canonical ensemble (i.e. the disk is 
allowed to exchange particles with a reservoir at infinity), in general there is no charge 
fluctuation and the charge Q is such that \Q\ has the smallest possible value (0 when 
possible, a fraction of the elementary charge q in the case of a one-component plasma 
with a background charge which is not an integer number of elementary charges — g). 
An exception is when the background charge of the one-component plasma is of the 
form— [(iV+ (l/2)]g with N integer. Then both values Q = ±(l/2)g are equally probable. 
We cannot explain why there is no charge fluctuation in the cases when the smallest 
possible value of the total charge is 0. Indeed, in these cases, bringing another elementary 
charge ±g from inflnity would cost only a flnite energy and one would expect Q = ±g to 
contribute to the grand canonical ensemble. That these values Q = ±g do not contribute 
might be a special feature of the solvable models at F = 2. We just do not know. 

Another tricky feature of two-dimensional Coulomb systems is that the two-dimensional 
Coulomb potential — ln(r/L) does not vanish at infinity, if the length L is finite. As dis- 
cussed in Appendix B, for obtaining sensible results in the Debye-Hiickel theory, it is 
necessary to take the limit L —>■ oo. 

APPENDIX A: DEBYE-HUCKEL THEORY IN A 
SPHERICAL OR CIRCULAR CAPACITOR 

In the spherical capacitor geometry, the solution of (|3.2| ) can be expanded in Legendre 
polynomials Pi{cos6) (where 6 is the angle between r and r') in the form 

oo 

KiT,T')=Y,ki{r,r')Pi{cos9) (A.l) 

1=0 

When this expansion is used in ( p.l|) and (|3.3|) , only the term / = survives the an- 
gular integrations. Thus, we only need the function ko{r,r'). The solution of (|3.2|) 
for an infinite system is K{r,r') = — /5exp(|r — r'|)/|r — r'| and its / = part is 
— /?sinh(Kr<) exp(— Kr>)/(Krr'), where r< (r>) is the smaller (the larger) of r and r'. 
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In the present geometry, there are additional "reflected" and "transmitted" terms. When 
the source point r' is in the inner sphere (r' < Ri), the solution is of the form 



ko{r, r') 



Krr 



y[sinh(fi;r<) exp(— Kr>) + asinh(Kr) sinh(/tr')] (r, r' < Ri] 



ko{r, r') = _^ ^i^h(KrO ^_&_ ^ (r' < R,, R, < r < R,) 



ko{r, r') 



/3d 



Krr 



- sinh(fi;r') exp(— Kr) (r' < Ri, r > R2) 



(A.2) 



This solution has the appropriate singularity at r = r', is otherwise regular at r = and 
r' = 0, and goes to when r — >■ 00. The four coefficients a, b, c, d are to be determined by 
the requirements that kQ{r,r') and dkQ{r,r')/dr be continuous at r = i?i and r = R2. In 
particular, one finds 



exp(— Ki?i)(/t_R2 — i^Ri 



[1 + KR2 — kRi) cosh(fi;/2i) + sinh(Ki?i 



(A.3) 



Using the first equation ( |A.2| ) and ( |A.3| ) in ( ^.11 ) and ( |3.3| ) gives 

In two dimensions, in the circular capacitor geometry, the calculation is very similar 
to the above one. The expansion ([A.lD is replaced by 



K{r,r') =J2kiir,r') cos{W) 
1=0 



(A.4) 



where the / = part is, in terms of modified Bessel functions Iq and Kq, when r' < Ri, 



ko{r, r') = -l3[IoiKr^)KoiKry) + a/o(/tr)/o(/tr')] (r, r' < Ri] 
ko{r,r') = -(3Io{Kr')[b\n{Kr) + c] (/ < Ri, Ri < r < R2) 
kQ{r, r') = -pdIo{Kr')KoiKr) (r' < R^, r > R2) 

By the same method as above, the coefficient a is found as 

~K,{kR{) + _ nRrK,{KR,) In f 



(A.5) 



(A.6) 



Using the first equation (|A.5|) and (|A.6|) in (|3.1|) and (|3.3|) , and the Wronskian of the 
Bessel functions, gives (|4.2|). 



APPENDIX B: DEBYE-HUCKEL THEORY IN A 
FINITE SYSTEM 

In Section 3 and Appendix A, the Debye-Hiickel differential equation ( |3.2| ) was writ- 
ten and solved in the spherical capacitor geometry; with minor modifications, the same 
approach holds in the two-dimensional case, for a circular capacitor, as discussed in Sec- 
tion 4.2 and Appendix A. In these geometries, the Coulomb fiuid extends to infinity in 
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the region r > R2 and it is obvious that the boundary condition should be K{y^ r') 
as r — s> 00) (In this approach perfect screening is globally satisfied: The charge of one 
particle plus the charge it induces in the two conductors sum to zero. However, there 
is no perfect screening if one takes into account only the charge of a particle sitting on 
one of the conductors and the charge it induces on that conductor only, and this gives 
rise to a charge fluctuation on each conductor). The case of one sphere or disk in the 
grand canonical ensemble was obtained by taking the limit R2 ^ 00 aX the end of the 
calculation. 

The question arises of how to formulate the Debye-Hiickel theory in the grand canonical 
ensemble, directly starting with only a sphere or a disk of radius R. What is the boundary 
condition to be imposed dX r = R 1 Choquard et al.p5|] have already investigated this 



problem. Nevertheless, we revisit it, hoping to clarify some delicate points. 

We start with the three-dimensional case. An unambiguous way of formulating the 
Debye-Hiickel theory is to start with a full diagrammatic expansion ||26|| in the grand canon- 



ical ensemble, and to make a topological reduction replacing the fugacity (fugacities) by 
the density (densities). The Debye-Hiickel correlation function is obtained by resumming 
a class of diagrams (the chain diagrams), or, equivalently, by taking for the function K 
in eq.(|3.lD the solution of the integral equation 

K{r, v') = --^ -ff dr"-^K{r", v') (B.l) 
|r — r I An J |r — r '| 

This integral equation can also be seen as the Ornstein-Zernicke equation in which the 
direct correlation function between particles of species a and b is approximated by —(3 
times their bare Coulomb interaction qaqb/\r — r'\. In a finite system, the densities in 
= AnP J2a ^a^ll are position-dependent near the boundaries. However, for the large-size 
systems considered here, this effect can be neglected and will be taken as a constant 
in the whole system. 

By taking the Laplacian of both sides of the integral equation (|B.1|) , one obtains the 
partial differential equation ( |3.2| ). However, the integral equation provides the boundary 
condition to be used in (|3.2| ). In the presently studied case of a sphere of radius R, we 
can use the expansion ([A.ID in Legendre polynomials. For brevity we only consider the 



I = part. The integral equation (p.l| ) gives for the boundary condition on the surface 
r = R (with r' < R) 

ko{R,r') = -^[1 + ^ rdr"r"'W,r')] (B.2) 
R p Jo 

Using in ( [B.2| ) the general form ( |A.2| ) of the solution of the partial differential equation 
( p.2|) (with now r, r' < R) and performing the integral determines the coefficient a as 

exp{-KR) 
cosh[KR) 

(The square bracket in the r.h.s. of ( [B.2|) does not vanish: there is no perfect screening 
on the sole sphere). This directly obtained result ( [B.3| ) is identical to the limit of the 
spherical capacitor a coefficient ( |A.3| ) (with Ri = R) as R2 —>■ 00. One recovers the mean 



14 



square charge (|2.4]) . The existence of a fluctuation confirms that the formulation of the 
Debye-Hiickel theory by the integral equation ( [B.l|) is a grand canonical one. 

The two-dimensional case is more tricky. The bare Coulomb interaction between two 
unit charges now is — ln(|r — r'|/L), with L some fixed length. The integral equation now 
is 

I 'I 2 I If \ 

K{v, r') = /3 In ^IZZl + — [ dr" In ~/ ^ K{r", r') (B.4) 
L 2tt J L 



with = ^TiPY.a^aQa- ^hc / = part of the expansion (|A.4|) obeys the boundary 
condition (with r' < R) 

koiR,r') = P\n^[l + ^ rdr"r"A;o(r",r')] (B.5) 
L p Jo 



Using in ( p.5| ) the general form ( |A.5| ) of the solution of the two-dimensional analog of 



the partial differential equation (|3.2| ) (with now r, r' < R) and performing the integral 
determines the coefficient a as 

kRKi{kR) \n^ + Ko{KR) 
" KRIi{KR)\n^ - Io{kR) ^ ■ ' 

This result ( p.6|) is unacceptable as it stands, since it depends on the arbitrary length 
L, which only determines the zero of the potential and should not enter the correlation 
functions. Actually, using a Coulomb interaction — ln(r/L) which does not vanish at 
infinity causes difficulties at many places. For instance, the Coulomb energy of a macro- 
scopic disk of radius R, carrying the macroscopic charge Q near its circumference, would 
be (1/2)Q^ ln(L/i?), negative if R> L. Thus, in the grand canonical ensemble, configu- 
rations of infinite \Q\ would be favored, causing the grand canonical partition function to 
diverge. This seems to indicate that the limit L ^ oo (such that the zero of the potential 
recedes to infinity) should be taken in ( [B.6| ), which then becomes 

This result ( [B.7| ) is identical to the limit of the circular capacitor a coefficient ( |A.6| ) (with 
Ri = R) as R2 ^ oo. It might even be noted that, if in ( A. 6 ) we take Ri = R and 



R2 = L ^ R and neglect the term of order R1/R2 = R/L in both the numerator and the 
denominator, ( p.6|) is recovered. This is a further indication that the limit L ^ 00 should 
be taken in the case of a system made of one disk only. Now, there is perfect screening 
(the square bracket in the r.h.s. of ( p.5|) vanishes) and there is no charge fluctuation on 
the disk. 

Still another way of dealing with a finite system, in three or two dimensions, would 
be to first assume that the whole space external to the system is filled with a medium 
of Debye wave number k' and solve the partial differential equation ( |3.2| ) (or its two- 
dimensional analog) in the whole space, taking /t^(r) = in the system and K^(r) = k'^ 
outside, with the proper continuity conditions on the system boundary and K{v, r') 
as r — > 00. One recovers the same results as above. 



15 



APPENDIX C: TWO-DIMENSIONAL 
TWO-COMPONENT PLASMA AT P = 2 



For this exactly solvable model[|l4[, the charge correlation function can be expressed 
in terms of Green functions (^^^(r, r') and G r') as 

(p(r)p(r')) = -2mV[|G++(r,r')P + |G_+(r,r')P] +n(r)g25(r-r') (C.l) 

where m is a rescaled fugacity such that (2m) ~^ is the bulk correlation length, and n(r) 
is the total number density (actually, n is a divergent quantity, but it will be seen that 
this divergence causes no difficulty here). In ( |C.1| ), the charge symmetry of the model 
has been taken into account. We consider the circular capacitor geometry and the source 
point r' is in the inner disk (r' < Ri). For (|C.1D to be the charge correlation function, r 
must be in a filled region (r < i?i or r > R2). 

When r is in a filled region, the Green functions obey the equations M 



(m^ - A)G++ (r, r') = m5{r - r') (C.2) 

and 

G_,(r,rO = -^f^ + i^)G,,(r,r') (^3) 

m \or r oip J 

where (r, 9?) are the polar coordinates of r. The solution of (|C.2|) is an expansion of the 
form 

m °° 

G++(r, r') = — V [Ii{mr<)Ki{mry) + aili{mr')li{mr)] exp[i/(v9 - ip')] (r, r' < Ri) 

/vr , 

(=— 00 

m °° 

G'++(r, r') = — E diIi{mr')Ki{mr) exp[i/(<^ - if')] {r' < R,, r > R^) (C.4) 
/vr , 

These expansions have the proper singularity at r = r', are otherwise regular at r = 
and r' = 0, and go to zero when r ^ 00. As to ( |G.3| ) gives 

G_+(r,rO = 

E [Ii{mr')Ki+i{mr) - aili{mr')li+i{mr)] exp[i(/ + l)ip - Hip'] {r < r < Ri) 

=—00 



2n , 



— diIi{mr')Ki+^{mr) exp[i(/ + 1)(^ - il^p'] (r' < r > R2) (C.5) 

ZTT , 

j=— 00 

When r is in the empty region Ri < r < R2 (where m = 0), as functions of r, 6*++ 

depends only on z = rexp{iip) and G ^ depends only on ^ = rexp(— iyj). Thus, the 

expansions are of the forms 

m °° 

G++(r,r') = — E biliimr'){mryexp[il{ip-ip')] {r' < R,, R, < r < R^) (C.6) 

(= — 00 
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and 



G_^(r, r') = — V Qj/(mr')(mr)-('+^) exp[i(/ + l)^ - il^'] (r' < R^, < r < R2) 

(C.7) 

The coefficients ai,bi,ci,di are to be determined by the requirements that 6*++ and 

G ^ be continuous at r = i?i and r = R2. In particular, after having used some functional 

relations between Bessel functions, one finds 



di 



lIi{X)Ki{aX){a^ - a'^) - XIi{X)K[{aX)a^+^ + XI[{X)Ki{aX)a 



-I 



-1 



(C.8) 



where X = mRi and a = i?2/-Ri- 

For computing the mean square charge on the inner disk, using the perfect screening 
relation / dr(p(r)p(r')) = 0, where the integral is on the whole space, it is convenient to 
rewrite (13.31) as 



dr 



dr'(p(r)p(r )) 



(C.9) 



'r>R2 Jr'<Ri 

and to use ( |C.1|) , omitting the self part nq'^6{r — r') which does not contribute to ( |C.9| ). 
The angular integrations are easily performed, by using the mutual orthogonality of the 
functions exp{ilip), and (p.9|) becomes 



(g2) = 2mV y / dr'r'l^{mr') / dr r[Kf (mr) + K^.^{mr)] (C.IO) 

where di is given by ( |C.8| ). After having performed the integrations in (|C.10| ) and used 
some functional relations between Bessel functions, one obtains ( ^.7] ). 



APPENDIX D: TWO-DIMENSIONAL 
ONE-COMPONENT PLASMA AT P = 2 



For this exactly solvable modeipD|,p2|, the charge density (including the background 



contribution) and the charge correlation function can be expressed in terms of a projector 
P(r,r') 

(p(r)) = g[P(r,r)-n] (D.l) 

and 

(p(r)p(r'))-(p(r))(p(r')) = g^[-|P(r,r')r + P(r,r)5(r-r')] (D.2) 

We consider the circular capacitor geometry. For (D.l) and (D.2) to be the charge density 
and correlation, respectively, r and r' must be in a ffiled region. 

The electric potential qV{r) created by the background will be needed. It obeys 
AV"(r) = 27m in the ffiled regions r < Ri and r > R2. It obeys AV{r) = in the empty 
region Pi < r < P2. At the boundaries r = Ri and r = P2, V{r) and d\^/dr must be 
continuous. Up to an overall irrelevant additive constant, 

V{r) = -Tcnr'^ (r < Pi) 
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V{r) = ^unRl + TinRl In ^ {Ri < r < R2) 
2 Ri 

V{r) = -Txn{Rl - Rj) + TxnRl + 
2 Ri 

l,rnr^ + l,rn(i??-i?2)ln^ {r > R2) 

L L 1X2 



(D.3) 



P(r, r') is the projector on the functional space spanned by the functions V^i(r) = 
exp[— V^(r)][rexp(i93)]^ (Z = 0,1,2,3,...) (this definition of "^i holds in the filled regions 
r < R\ and r > R2, while ipi{r) = in the empty region Ri < r < R2): 

^(r, r') = E 77^/(r)^z(r') = E 77 exp[-r(r) - V(r')]rV' exp[i/(v9 - ^')] (D-4) 
i=o 1=0 

where C; is the normalization constant 



Ci = {[ + f )dr exp[-2F(r)]r2^ 

Jr<Ri Jr>R2 

= zr4rT[7(^ + l,>^i) + ^(^ + l,n,r,)] 



(D.5) 



where the functions 7 and D are defined in ( [4.11 ) and ( 4.14 ). 
The average charge on the inner disk is 



(Q) = / dr(p(r)) 

Jr<Ri 



(D.6) 



Using ( |D.1| ) and ( |D.4| ) in ( p.6| ) and performing the integration gives ( |4.13| ) 
The charge fluctuation on the inner disk is 



dr / dr'[(p(r)p(r'))-(p(r))(p(r'))] 

r<Ri Jr'<Ri 



(D.7) 



Using ( p.2| ) and ( p.4| ) in ( p.7| ) and performing the integrations, using the mutual orthog- 
onality of the functions exp{ilip), gives ([4.15|) . 
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